Sound propagation in a Fermi gas near a Feshbach resonance 
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Sound waves are observed and studied in an optically trapped degenerate Fermi gas of spin-up 
and spin-down atoms with magnetically tunable interactions. Measurements are made throughout 
the crossover region, from a weakly-interacting Fermi gas through the resonant Fermi superfluid 
regime to a Bose condensate of dimer molecules. The measured sound velocities test the equation 
of state and confirm the universal hypothesis. 



Optically-trapped, strongly-interacting Fermi gases [l| 
provide a unique laboratory for testing nonperturbative 
many-body theories in a variety of fields, from neutron 
stars and nuclear matter Q to quark-gluon plasmas Q 
and high temperature superconductors In contrast 
to these other systems, interactions in ultra cold atomic 
Fermi gases are tunable using the Feshbach resonance 
phenomenon At magnetic fields well above the reso- 
nance, one obtains a weakly attractive Fermi gas, while 
well below the resonance, spin-up and spin-down atoms 
are joined into dimer molecules, which form a Bose- 
Einstein condensate (BEC) 0, 0, 0, ij. On resonance, 
the gas is a strongly interacting Fermi superfluid [Io| . 

The focus of this Letter is the precision measurement 
of the sound velocity at low temperatures in this unique 
quantum system. The experiment addresses the rich 
physics of wave propagation in a strongly-interacting sys- 
tem, such as possible coupling between first and second 
sound 11] and nonlinear behavior. We find that the 
wavefront is flat despite nonuniform density across the 
wavefront. This observation verifies the predicted rela- 



tion 12| between the measured velocity in a cigar-shaped 
trap and the equation of state. As a result, by mea- 
suring the velocity we probe zero-temperature equations 
of state niacin), where ni oc is the local chemical poten- 
tial and n the total density of atoms. The equation of 
state is a central result of nonperturbative many-body 
theories that determine the thermodynamic parameters. 
Measurements at resonance over a wide range of densities 
confirm the universal hypothesis [j], 0, [Til. 14 1 as well as 
the location of the Feshbach resonance [l5j . 

The medium for the sound propagation is a 50:50 mix- 
ture of 6 Li atoms in the two lowest internal states or a 
BEC of Li2 molecules near a broad Feshbach resonance 
centered at B = 834 G [HI]. The gas is cooled by forced 
evaporation in an ultrastable CO2 laser trap at the res- 
onant magnetic field [l[ • Near the end of evaporation, B 
is adiabatically set to a chosen final value between 700 
and 1100 G, and the trap is adiabatically recompressed 
to a desired depth Uq between 0.14 and 12 fiK. 

The experimental samples are nearly in the ground 
state, which can be determined from the absorptive im- 
ages of the gas [ij. The gas is imaged after release from 
the trap and expanded to avoid saturation of images 



due to high absorption. Close to unitarity and on the 
Fermi side of the Feshbach resonance, the density pro- 
files can be closely approximated by the zero-temperature 
Thomas-Fermi profiles of a non-interacting Fermi gas. In 
the molecular BEC regime, the clouds do not have any 
observable thermal component, which appears when the 
cooling efficiency is significantly reduced. From the ab- 
sorptive images we measure the total number of atoms 
N in the range 0.6 to 4.9 x 10 5 . 

A sound wave is excited in the sample by the same 
method as used for a BEC of atoms [lgj. A thin slice 
of green light bisects the cigar shaped cloud. The green 
light at 532 nm is blue detuned from the 671 nm transi- 
tion in lithium, creating a knife which locally repels the 
atoms. The 1/e thickness of the knife is ~ 30 /im. The 
height of the knife potential, Uknife, is estimated from 
the size, power, and atomic polarizability. Uknife is cho- 
sen to be between 0.5 and 14 fiK. The knife is pulsed on 
for 280 /xs, much shorter than typical sound propagation 
times ~ 10 ms. The resultant perturbation propagates 
outward along the axial coordinate z. We allow propa- 
gation for a variable amount of time and then release the 
cloud, let it expand for a fixed time, and image destruc- 
tively. Propagation of the density perturbation is shown 
in Figure Q] The propagating feature consists of a region 
of low density (valley) with regions of high density (peak) 
on either side. 

The main result of the Letter, shown in FigureEJ is the 
measurement of co/vp, the normalized sound velocity in 
the center plane (at z — 0) vs 1/kpa, the dimensionless 
interaction parameter. The reference Fermi velocity vp 
is that of a non-interacting Fermi gas at the trap cen- 
ter, i.e., E F = hiuJxLUyUJz) 1 / 3 (3A) 1 / 3 = m atom v 2 F /2 = 
h 2 k F /2m a tom, where Ep is the Fermi energy and kp 
is the corresponding Fermi wave vector. The trap fre- 
quencies tOi are individually calibrated for each velocity 
measurement. At zero temperature, 1/kpa is the only 
quantity that determines the physics. The s-wave scat- 



tering length a — a(B) is estimated from Ref. [la]. At 
\kFd\ 1 the system is in the strongly interacting regime 
and is a superfluid. The BEC (weakly interacting Fermi) 
regime corresponds to 1/kpa > 1 (1/kpa < —1). We 
tune kpa by changing B and Uq. 

Five measures are taken to reduce systematic errors. 
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FIG. 1: Sound propagation at the Feshbach resonance: Inte- 
grated axial profiles of the cloud for different in-trap propaga- 
tion times. Positions of the propagating density depressions 
(valleys) are highlighted. Propagation times are shown. 
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FIG. 3: Axial positions z of the left- and right-traveling den- 
sity features of [T]/s the in-trap propagation time. Coordi- 
nates of the valleys (blue-closed) and front peaks (red-open) 
are shown. Dashed lines correspond to the measured Fermi 
radii (z = ±R Z ). 



z. Making a general approximation for the equation 
of state fJ,i oc (n) — C n 1 and calculating the velocity 
within using hydrodynamic theory [13 , Il9| . one finds 
that c(z) oc y/l — z 2 /R 2 for a harmonically trapped 
gas, where R z is the axial Thomas-Fermi radius. Then, 
one may show that a small density perturbation prop- 
agates as z(t) = R z sm(±cot/R z + if), where the phase 
ip = arcsm(zknife/ Rz) depends on the excitation posi- 
tion. Co, the velocity of the features at z = 0, is obtained 
by fitting the coordinates of the propagating features 
which lie within the central 65% of the cloud (Fig.[3|). 

Fourth, we account for nonlincarity in propagation: 
The front peak moves faster than the valley, as predicted 
earlier for BECs • In Fig- HI we show the velocities of 
the valley and peak vs the excitation strength, Uknife/ ' Hi 
where fi is the global chemical potential estimated us- 
ing Monte-Carlo modelling 171 l2l| ■ The peak and valley 
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FIG. 2: Sound velocity in the center plane normalized to vp 
vs the interaction parameter, l/fej?a. Black dotted curve - 
mean-field theory that uses the Leggett ground sta te I 17ll . Red 
solid curve - quantum Monte Carlo calculation [17T |. Black 
dashed curve - Thomas-Fermi theory of a molecular BEC, 
using a mo i = 0.6 a. Open and closed circles represent ranges 
of trap depths, 0.6 to 12 fjK and 140 to 410 nK, respectively. 



First, by tracking both right- and left-propagating fea- 
tures (Fig. [3]) we eliminate effects due to the overall mo- 
tion of the cloud. Second, we correct for a small change 
of the axial cloud size between the release and imaging 
(< 4.7%, 1.2% on average). In order to obtain in-trap 
positions of the propagating features, each imaged po- 
sition is scaled by a coefficient of expansion calculated 
within a hydrodynamic model [l], [l8[ . 

Third, we account for the dependance of the ve- 
locity on the non-uniform density distribution along 
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FIG. 4: Valley (red-closed) peak (blue-closed) and mean 
(green-open) velocity vs excitation potential normalized to 
the global chemical potential. Measured at the Feshbach res- 
onance. 

velocities vary with excitation strength, but the mean 
velocity remains constant. All velocities converge to a 
single value at low excitation strength showing that the 
mean accurately represents the velocity of an infinitely 
small perturbation. In reporting data, we use the mean 
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value. For data of Fig. [2l Uknife/^ ranges from 1 to 30. 
We do not observe formation of shock-waves which are 
predicted for a BEC [U. 

Fifth, the data is corrected for anharmonicity of the 
trapping potential. We correct the measured sound ve- 
locities and find the velocity values which would be mea- 
sured in a perfectly harmonic trap. Only the radial an- 
harmonicity is taken into account because for shallow 
traps, in which the correction is greatest, harmonic mag- 
netic confinement dominates axially. The radial confine- 
ment is optical, and the potential is gaussian. Axially, the 
potential is a sum of an optical and a parabolic magnetic 
potential with a common center (oj^ = uJ^ opt + B ). 
The corrections for cq have been calculated to the first 
order in n/Uo from Eq. [T] below using either a Fermi 
(7 = 2/3) or Bosc (7 = 1) equation of state fii oc cx rO . 
The respective corrections are co/c™ eas ' = 1 + 0.10 h/Uq 
and c /c^ eas - = 1 + 0.14 fi/U . The Fermi and Bose cor- 
rections have been used for l/kpa < and l/kpa > 1. 
For < l/kpa < 1, these two corrections were linearly 
combined. pi/Ua ranges from 0.025 to 0.43. 

Precision characterization of the trap is important for 
calculating the reference Fermi speed vp. The mean ra- 
dial frequency ui± = ^uj x uj v is calculated from the radial 
breathing mode frequency lObreath j22, 23] , which is mea- 
sured for each point of Fig. [2] on the day of the respec- 
tive sound measurement. LObreath is measured either at 
resonance, where oJbreath = \/10/3oi;j_, or in the far BEC 
regime, which yields consistent results. In the latter case, 
we assume u>breath. =_2w 1 neglecting a small beyond- 
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The ratio u> x /ui y = 1.056 



and the total axial frequency lu z has been measured us- 
ing parametric resonance in a weakly-interacting Fermi 
gas 22[. The frequency of the magnetic field poten- 
tial, ujb/2tt = 20.5 Hz at 834 G, is measured using a 
sloshing mode. The trap anisotropy ratio lu±/ui z is be- 
tween 4.1 and 25 depending on Uq. For Uq = 6.1 /iK, 
w ± /2tt = 690(2) Hz and u z /2-k = 34.2(0.2) Hz. All 
radial frequencies have been corrected for anharmonic- 



ity 



23l l26l|. The individual corrections to cq and vf 
partially cancel when calculating the ratio co/vp. The 
resulting correction to cq/vf is within 3.1% and, on av- 
erage, is 0.7%. 

Linking the speed of sound to the equation of state 
enables quantitative tests of theories. By analyzing the 
wave front shape we now show that the hydrodynamic 
model of Capuzzi et al. [12], for a cigar-shaped trap, pro- 
vides a correct description of propagation. The primary 
result of this model is that the wave front stays flat, as 
described by 
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ITT'atom J (<9/i; oc /<9n) dx dy 



(1) 



In contrast, a simple model of isotropic sound propaga- 
tion assumes ci oc (n) — \J (dP/dn) / m a tom, where P is the 



pressure. In this case, the wave front curves as the wave 
propagates. For the same equation of state, these two 
theories predict a few per cent difference in the speed of 
the perturbation. In Figure [5] we experimentally demon- 
strate that the speed is independent of the transverse 
position, which validates Eq.[T]for our experimental con- 
ditions. And, hence, we can rely on Eq.[T]in connecting cq 
to (J,ioc(n). At this point we are in a position to compare 
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FIG. 5: Local sound velocity vs transverse coordinate x, in 
units of the Thomas-Fermi Radius, at the feshbach resonance 
showing flat wavefront. Images were divided into 21 longitudi- 
nal segments of equal width. For each ith segment, a(z = 0) 
has been measured. The dashed curve represents the line of 
site averaged isotropic sound speed. 

the data of Fig. [2] to predictions. 

In the BEC regime, at l/kpa > 1, the experiment tests 
the interaction properties of composite bosons, molecules 
made out of two fermions. One may approximate the 
state of the g BEC of singlet molecules that col- 

lide with s-wave scattering length a mo /. Using the lo- 
cal equation of state (J>ioc(n mo i) = 2nh 2 a mo in/ m mo i in 
Eq. [U one finds co/vp — (5/2 /ci?a mo /) 1 / 5 /4. An exten- 
sion of the Leggett ground state into the Bose regime [27 1 
gives a mo i — 2a{B) resulting in a prediction for the 
speed 25 — 40% higher than measured. An exact so- 
lution of the four-fcrmion scattering problem predicts 
dmoi — 0.6 a [ll] yielding a speed (black dashed curve 
in Figure [2]) close to the measurements. 

Theories covering the whole crossover from the Fermi 
to Bose regime can be tested against our data. In Fig. [21 



the black dotted curve [17 [represents the theory based on 
the Leggett ground state [2jJ . The theory predicts signif- 
icantly higher speed than measured except for the Fermi 
side of the resonance. The red solid curve [13] shows the 
calculation that uses the equation of state obtained from 
a Monte-Carlo simulation 2l| . This theory shows much 
better agreement with the data. Interestingly, in most 
of the molecular BEC regime, at l/kpa > 1, the data is 
systematically lower than the prediction. 

The lower measured speed at 1 < l/kpa < 5 could 
be due to the coupling of first and second sound, which 
travels slower (ll| . High temperature cannot be a reason 
for lower speed. We have found that higher tempera- 
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tures correspond to higher sound velocities. It is also 
unlikely that the speed is lower due to a dark soliton, 
which travel slower than sound, and, in principle may 
exist in our molecular BECs. The tracked density depres- 
sions are most likely not dark solitons for three reasons. 
First, unlike a dark soliton, the feature spreads out and 
loses contrast with time (Fig.Q}, which is observed at all 
1/kpa's. Second, the soliton size is predicted to be of the 
same order as the coherence length [29j , h / y/2fi m mo i < 1 
/im, while the imprinted density depression has the width 
of 10 — 20 jum. Third, the optimum soliton excitation re- 
quires phase imprint of it [30(, while in our experiment, 
the phase imprint for most of the BEC-side experiments 
is ~ 20tt. 

The small speed increase at 1/kpa > 5 could occur due 
to a technical reason: To achieve these high 1/kpa values 
we have to use very shallow traps. The BECs released 
from these weakest traps are small and difficult to image. 

The universal hypothesis is positively tested in our 
measurements. The hypothesis states that at resonance, 
where kpa — oo, cq/vf should be independent of density. 
Indeed, within the experimental precision, at resonance, 
at 834(2)G, we observe that Cq/vf remains constant as 
the trap depth is changed from 410 nK to 12 /iK. In 
contrast, below the Feshbach resonance, at 821 G, co/vf 
decreases with decreasing Uq. These measurements of 
Cq/vf vs Uo show that the Feshbach resonance location 
is much closer to 834 G than to 822 G [Hj]. 

From the sound velocity at resonance, we have also 
measured the universal constant /? [H, [2] , which connects 
the local interaction and Fermi energy in the universal 
regime: U int = /3e F - Using /i ioc = e F + U mt = (1 + 
0) (37r 2 n) 3 / 2 h 2 /2m atom and Eq. [1] one may relate cq to 
(3 as: 
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from which we find /? = -0.570(0.015). 

On the far Fermi side of the Feshbach resonance, at 
1/kpa < —1.3, propagation of sound is not observed. 
Instead, the hole from the laser knife fills without propa- 
gating, as in a non- interacting gas. At the temperatures 
we achieve, the gas is not necessarily a supcrfluid in the 
weakly interacting Fermi gas region, and hydrodynamic 
sound propagation is provided by collisional, rather than 
superfluid hydrodynamics. 

In conclusion, we have observed sound propagation, 
connected observations to a hydrodynamic model, and 
provided data for testing theories describing strongly- 
interacting Fermi systems as well as their crossover into 
the weakly interacting Fermi and Bose regimes. 
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